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We determine tlie structure of cliarge-stabilized colloidal suspensions at low ionic strength over 
an extended range of particle volume fractions using a combination of light and small angle neutron 
scattering experiments. The variation of the structure factor with concentration is analyzed within 
a one-component model of a colloidal suspension. We show that the observed structural behavior 
corresponds to a non-monotonic density dependence of the colloid effective charge and the mean 
interparticle interaction energy. Our findings are corroborated by similar observations from primitive 
model computer simulations of salt-free colloidal suspensions. 
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Deionized colloidal dispersions have recently attracted 
much experimental [J, 0, Q and theoretical attention 
d, [^, [1, 0] because of their unique role as simple model 
systems in the field of surface charge regulation. In addi- 
tion to experimental challenges posed by these systems, 
such as the need of control over the ionic strength of 
the medium, a number of intriguing observations have 
been made that challenge classical theories of electro- 
static screening and electrokinetics (e.g., see and 
references therein) . Some of the conceptual difficulties in 
the understanding of these systems stem from the coun- 
terion dominance in the ionic cloud of the colloid, which 
makes the screening ability of the medium dependent on 
the colloidal concentration. At the same time, the coun- 
terion condensation, which is controlled by the balance 
of the counterion chemical potential between the bulk 
suspension and the condensed layer, would depend on 
the concentration through the state-dependent screen- 
ing parameters [1, [13, EH- Therefore, modeling these 
systems within the usual one-component description can 
only be achieved using effective density-dependent inter- 
action parameters, charge and screening length . The 
dependence of the colloid effective charge on the volume 
fraction 0, [^, or background ionic strength Q has 
been predicted in a few mean-field models and is expected 
to be non-monotonic. However, this non-monotonic be- 
havior has never been validated experimentally. 

In this Letter, we present a study of the density de- 
pendence of the colloidal structure and interactions at 
low ionic strength. Experiments have been carried out 
over a wider range of concentrations than in any pre- 
vious work. We analyze the particle structure factor 
obtained from light and neutron scattering at different 
particle volume fractions and demonstrate that (i) the 
variation of colloidal structure in deionized suspensions 
can be interpreted in terms of a one-component model 



(OCM) that considers Yukawa-like pair interactions be- 
tween colloids with state-dependent parameters and, (ii) 
the non-monotonic concentration dependence of the in- 
teraction energy, as characterized by the colloid effective 
charge, is a general feature of deionized colloidal disper- 
sions. Primitive model (PM) numerical simulations of 
salt-free suspensions corroborate our observations. 

Stock suspensions of negatively charged particles (sul- 
fonate polystyrene latex) were purchased from Interfacial 
Dynamics Corporation (Portland, USA) and master dis- 
persions were prepared following the protocol described 
in Ref. 



mil 



The samples were prepared by dilut- 
ing the master suspensions directly into quartz cells with 
water/ethanol mixtures; filtered several times (pore size 
0.2 /im) to remove dust particles. These mixtures were 
chosen in order to avoid crystallization by modifying the 
solvent dielectric constant [ITj] , thereby providing a well 
controlled model system of charged dispersions with vol- 
ume fractions as high as ^ 16%. Accurate values of parti- 
cle radii and polydispersities, as obtained from DLS and 
SANS experiments, are given in Table HI 

The effective interparticle interaction can be quanti- 
fied through the analysis of the suspension static struc- 
ture. This requires probing length scales comparable 
to the size of the particles. Usually such information 
can be obtained from static and dynamic light scatter- 



TABLE I: Polystyrene latex spheres used in the experiments. 
Bjerrum length, Ab = / (47r e ep fcs T) , from the dielectric 
permittivity data in Ref. [13] • 



Sample # 


Radius 


Bjerrum Volume Scattering 




(nm) 


length (nm) fraction (%) Technique 


SI 


54.9 ± 6.4 


1.296 0.035 - 01.15 3D-DLS 


S2 


58.7 ± 8.6 


1.383 0.135 - 01.48 3D-DLS 


S3 


54.7 ± 5.3 


1.482 2.500 - 15.50 SANS 
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ing (SLS/DLS) in the single scattering regime [l^|. In 
dense suspensions, however, these techniques cannot be 
apphed in most cases due to multiple scattering of light. 
An elegant way to overcome this limitation is the use of 
modern scattering techniques such as 3D dynamic light 
scattering (3D-DLS) [l^ or small angle neutron scatter- 
ing (SANS). While 3D-DLS provides valuable informa- 
tion at small and intermediate values of the scattering 
vector q, access to large g-vectors from SANS allows to 
normalize the static structure factor S{q) and, in that 
way, unambiguously determine the suspension structure. 
A detailed description of the light scattering (LS) exper- 
iments is given in Ref . [l5| . 

Multiple scattering of neutrons in our SANS experi- 
ments has been suppressed by partially contrast match- 
ing the particles using H20/ethanol - D2 0/deuterated 
ethanol mixtures. The samples were kept in stoppered 
quartz cells (Hellma, Germany) with a path length 
of 2 mm and containing a mixed-bed of ion-exchanger 
resins; deionization was typically completed within 2 
weeks. The measurements were performed with a mean 
neutron wavelength of 1.27 nm and at a detector dis- 
tance of 20.3 m, which corresponds to a g-range of 0.01- 
0.1 nm^^. For details on the initial data treatment and 
data analysis see Ref. [l^l and references therein. The 
SANS scattering intensity can be expressed as I{q) — 
AP(g)sANS '5'(g)sANS, where P(g)sANS is the normal- 
ized effective particle form factor, S'(g)sANS is the effec- 
tive interparticle structure factor and A is an amplitude, 
which is proportional to the particle number density, n, 
and to the square of the neutron scattering length con- 
trast between the particles and background fluid 21 1. 
The effective SANS static structure factor can be de- 
termined from S'(q)sANS = {Ao/A)[I{q)/Io{q)], where 
Io{q) = Aq P{q)sANS is the intensity scattered by a sus- 
pension of non-interacting particles with amplitude Aq. 

On the theory side, the static structure factor 
of colloidal dispersions can be calculated by solving 
the Ornstein-Zernike (OZ) equation together with the 
Rogers- Young (RY) closure relation fi^ for particles in- 
teracting with an effective pair potential Ues{r) of the 
Yukawa form 



/3ucff (r) = Z^ff Ab 



exp(Koff' a) 



1 + Kcff a 



exp(-Kcffr) 



(1) 



with an effective charge Z^g, colloid radius a, ther- 
mal energy = k^T, screening parameter = 
47r Ab {Zcsn + Us) and concentration of salt ions Us = 
2 X lOOOiVACs, where Cs is the molar concentration [l3 |. 
This procedure has been successfully used for calculating 
the structure in the strong screening regime at KcS a » 1 
[2^ but it becomes less accurate in systems with thick 
double layers as the interparticle potential deviates from 
the Yukawa shape [23[ . Despite this deficiency, an appro- 
priate choice of ZoS still allows one to predict the struc- 
ture factor of strongly interacting systems over a wide 



range of conditions as long as the pair interaction energy 
at the mean interparticle distance is correctly reproduced 
[i^ l. Nevertheless, an agreement of the OCM structure 
and the mean interaction energy with the original ones 
does not guarantee that other thermodynamic proper- 
ties, such as the osmotic pressure, can be predicted using 
the effective parameters. The difficulty with the pressure 
originates from the loss of the ionic degrees of freedom in 
the coarse-grained description. Various solutions to this 
issue have been suggested recently. In one, the missing 
pressure contributions have been recovered via inclusion 
of the boundary terms Q. Alternatively, one can use 
a modified OZ-RY approach to get the effective parame- 
ters in the fit procedure |5|]. The latter replaces the OCM 
(Yukawa fluid) pressure by the osmotic pressure known 
a priori from primitive model simulations or mean-fleld 
models (e.g. jellium) in calculating the isothermal com- 
pressibility. Wc use this route to extract the effective 
charges and effective potentials from the structure data 
measured in experiment and simulations. 

Fig. [1] displays the normalized scattering intensity for 
two of the samples S3 in Table [J a) cp = 2.5% and b) 
(p = 7.6%. In a), the squares are from 3D-DLS and 
circles from SANS measurements, thus illustrating the 
complementarity of both experimental techniques. The 
experimental error is estimated from the intensity statis- 
tics in a series of measurements. For both SANS and LS, 
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FIG. 1: (color online) Normalized scattering intensities of 
samples a.) (p — 0.025 and, h) (p — 0.076. The squares are 
from 3D-DLS and the circles from SANS. The dashed line is 
the calculated S{q) and the continuous line its SANS exper- 
imental convolution, 5'(g)sANS. Insets: scattering intensities 
for a) a diluted screened suspension and, h) p = 0.076. 



3 



the experimental data are in good agreement with the- 
oretical calculations (lines) based on integral equations 
theory using the polydisperse RY closure relation [2§| 
as described in Ref. The background salt concen- 
tration, Cs — 2 X 10~^M, was determined by measuring 
the conductivity of our deionized solvent mixtures as in 
Ref. 25|. Thus, the free parameters in our fitting proce- 
dure are the volume fraction if and the effective charge 
Zctf- The dashed lines represent the calculated S{q) and 
the continuous lines correspond to the result after its ex- 
perimental convolution, S'(g)sANS; with the appropriate 
resolution function [i^. The structure factor exhibits a 
pronounced maximum, whose positio n qm ax varies with 
the particle volume fraction as n"^/"^ [16|. Additionally, 
in the Insets we also show the SANS scattering intensi- 
ties for: a) a dilute, f ~ 0.005, screened (cg ~ 0.005 M 
of KCl) suspension and, b) the sample (p = 7.6%. Here, 
the lines are experimentally convoluted fits with poly- 
disperse form, Io{q) oc P{q)sA'NS, and structure factor, 
I{q) oc P(g)sANS 'S'(g)sANs- 

In Fig. [2^) we present the (/3-dependence of the height 
of the main peak of the structure factor, S'(qmax), as ob- 
tained from the LS and SANS analysis. At low volume 
fractions, the peak height is about 2 and varies weakly 
with concentration until ip « 0.002, then it decreases in 
the range 0.002 < < 0.01 to finally grow again for 
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FIG. 2: (color online) a) S'(qmax) vs. p. Squares and cir- 
cles are from 3D-DLS measurements in SI and S2, respec- 
tively. Triangles are from 3D-DLS (solid) and SANS analysis 
(open) in S3. Inset: simulations from PM (solid circles) and 
OCM-NPT (open circles), b) The colloid interaction energies 
/3«cff(»') at the mean interparticle distance Vm for experiments 
and simulations (Inset). Lines are guides to the eye. 



(fi > 0.01. We note that for > 0.06, the values of 
•^(gmax) > 3.0 indicate a glassy state [27[. While the 
increase of S'(gniax) with concentration can be readily ex- 
plained by shortening the distance between neighboring 
colloids and rising average interparticle repulsion, the de- 
crease at « 0.01 is a more complex feature. This de- 
pression can be attributed to the maximum of ionic con- 
densation or a minimum of the effective macroion charge 
as discussed below. To test whether this behavior is 
specific to the experimental systems considered here, in 
the Inset of Fig. [2^) we show results of PM simulations 
(solid symbols) of a salt-free charged colloidal suspen- 
sion with charge asymmetry 200 : 1 . The simulations were 
performed for systems with Ae/a = 0.2 using the same 
cluster Monte Carlo simulation protocol and settings as 
in Ref. ^2§\ with 80 colloidal particles. The 5'(ginax) from 
simulations shows a similar trend to the experiments al- 
though the main features are less pronounced. Never- 
theless, one can show that this characteristic trend of 
the structure factor corresponds to a quite peculiar be- 
havior of the pair interaction potential. Fig. [^b) de- 
picts the interaction energy between nearest neighbor 
colloids, Puef{{r^), as extracted from fits to the struc- 
ture data of systems S1-S3. /3ucff(7'm) follows the trend 
observed for S'(gmax) with more pronounced features: a 
maximum at </? « 0.001, a minimum at « 0.01, and 
an increase at higher volume fractions. The same holds 
for the PM results; solid circles in the Inset of Fig. [21d). 
Furthermore, we were able to self-consistently prove our 
primitive model results by using their associated effective 
charges and screening lengths in OCM-NPT simulations 
(open squares in the Inset of Fig. H^)). Here, we used 
N = 1000 colloidal particles and the osmotic pressure as 
obtained from the PM simulations. 

The interactions between colloids are conveniently 
characterized by their effective charge. To compare the 
results between different systems we use the dimension- 
less ratio ZcS — Z^S'^b /o. instead of the valency itself. In 
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FIG. 3: (color online) Z^s vs. ip from experiments (main) and 
simulations (Inset). The lines are calculations based on renor- 
malized jellium (dashed) and PB-cell (dash-dotted) models, 
and considering Cb — 2 x 10~^ M (main) and Cb — (Inset). 
Same symbols as in Fig. [2] 
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Fig. [3] we show the resuhs of ZcS as obtained from fitting 
the scattering data with our OZ-RY scheme (same sym- 
bols as in Fig. [2). Furthermore, in Fig. [3] we compare the 
data with results of the PM simulations described above 
(solid circles in the Inset), where ZbarcAe/a = 40.0 and 
the Zcfi are extracted following the same protocol as in 
the experiments. Both sets of data show a pronounced 
minimum of the effective charge at ip ~ 10~^. 

The behavior of Zeg in the range cp < 10^^ stems 
from the variation of the mean counterion concentration. 
Upon increasing ip, counterions in average get closer to 
the colloidal surface thus leading to a gradual decrease 
of ZcB- At (f > 10~^, in contrast, the effective charge 
starts increasing because the average electrostatic poten- 
tial in the bulk, or at least at r > r,ii, becomes com- 
parable to that on the colloid surface thus making the 
counterion condensation less favorable. In Fig.[3l we also 
present results from two models, which predict a similar 
density dependence for ZcS'. renormalized jellium Q and 
Poisson-Boltzmann cell model (PB-cell) We used 

Cs = 2 X 10^^ M in the experiments (main figure) and 
Cs = in the simulations (Inset). These models cor- 
rectly reproduce the qualitative behavior of the effective 
charge, although the values corresponding to the satura- 
tion regime are too high. This quantitative disagreement 
might be caused by various reasons. First, the best fit 
in our procedure corresponds to slightly lower effective 
charges than those in the jellium model. Moreover, the 
values are affected by the actual surface charge density 
on the colloids in experiment and ionic correlation effects 
in the simulations [28| . The charge dissociation in water- 
ethanol mixtures is expected to be lower than in pure 
water. Our present data suggest the dissociation on the 
level of ca. 20% of that in water, which means the actual 
charges are below the saturation values for these parti- 
cles at some volume fractions. Finally, we should note 
that the non-monotonic behavior of ZcS is characteris- 
tic for the countcrion-dominated screening regime. Our 
numerical estimations show that already amounts of salt 
comparable to the effective counterion concentration, ca. 

in our systems, might lead to disappearance of the 
minimum in the effective interaction and effective charge 
as function of colloidal volume fraction. 

In summary, we have studied the variation of the col- 
loidal structure factor in charged deionized colloidal dis- 
persions and demonstrated that its density dependence 
corresponds to a non-monotonic variation of the mean in- 
terparticle interaction and particle effective charge, with 
a minimum at about a volume fraction of 1%. Our re- 
sults cover a wider range of particle concentrations than 
in any previous study thus giving an excellent benchmark 
for models of electrostatic screening in colloidal disper- 
sions. Furthermore, our findings are confirmed by com- 
puter simulations in the counterion dominated regime of 
both primitive and one-component models. 
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